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Abstract 


The  Blaschke-Lebesgue  theorem  states  that  of  all  plane  sets  of 
given  constant  width  the  Reuleaux  triangle  has  least  area.  The  area 
to  be  minimized  is  a  functional  involving  the  support  function  and 
the  radius  of  curvature  of  the  set.  The  support  function  satisfies  a 
second  order  ordinary  differential  equation  where  the  radius  of 
curvature  is  the  control  parameter.  The  radius  of  curvature  of  a 
plane  set  of  constant  width  is  non-negative  and  bounded  above.  Thus 
we  can  formulate  and  analyze  the  Blaschke-Lebesgue  theorem  as  an 
optimal  control  problem. 


Introduct i on 


The  width  of  a  closed  convex  curve  in  a  given  direction  is  the 
distance  between  two  parallel  supporting  lines  perpendicular  to  that 
direction.  A  set  of  constant  width  b  has  the  same  width  in  all 
directions.  Besides  circle,  the  best  known  closed  convex  curve  of 
constant  width  b  is  the  Reuleaux  triangle  of  width  b,  i.e.  a  set  in 
R2  whose  boundary  consists  of  three  congruent  circular  arcs  of  radius 
b.  See  Figure  1.  The  Blaschke-Lebesgue  theorem  states  that  the 


Reuleaux  triangle  has  the  least  area  of  all  plane  convex  sets  of  the 


same  constant  width  b.  The  minimum  area  is  ^ T  ^  .  This  theorem 


was  first  proved  independently  by  Blaschke  [2] ,  and  Lebesgue  [17] . 
Besicovitch  [1] ,  Chakerian  [9] ,  and  Eggleston  [10; 11]  contain  a  proof 
of  the  Blaschke-Lebesgue  theorem. 

Optimal  control  theory  can  be  applied  to  geometric  extremum 
problems  for  plane  curves  as  follows:  The  functional  for  which 
extrema  are  examined  are  geometric  invariants  such  as  area  or 
perimeter.  The  system  of  ordinary  differential  equations  for  the 
control  theory  formulation  is  derived  from  the  Frenet-Serret 
formulas,  and  the  control  parameter  is  curvature. 

Klotzler  [16]  has  used  optimal  control  theory  to  study  n- 
orbiforms.  These  are  convex  planar  domains  which  can  be  rotated 
inside  a  regular  n-gon  under  tangential  contact  on  all  sides.  Plane 


sets  of  constant  width  can  be  rotated  inside  a  square  with  tangential 
contact  on  all  sides.  Our  approach  to  plane  sets  of  constant  w-idth 


i 

is  different  than  that  in  Klotzler  in  the  choice  of  the  functional  to 
be  minimized. 

In  the  following  we  discuss  preliminary  definitions  related  to 
sets  of  constant  width  and  include  necessary  background  from  the 
theory  of  optimal  control.  We  then  formulate  and  analyze  the 
Blaschke-Lebesgue  theorem  as  an  optimal  control  problem. 

I 

Prel iminaries 

By  a  convex  body  in  Rn  we  mean  a  compact  convex  subset  of  Rn  with 

i 

nonempty  interior.  For  each  direction  u  £  S'1’1,  where  S'1"1  is  the  unit 
sphere  centered  at  the  origin  in  Rn ,  we  let  h(K,u)  denote  the  support 
function  of  the  convex  body  K  evaluated  at  u.  Thus, 

I 

(1)  h(K,u)  =  sup{u-x:  x  £  K}  , 

which  may  be  interpreted  as  the  distance  from  the  origin  to  the 

: 

supporting  hyperplane  of  K  having  outward-pointing  normal  u.  The 
width  of  K  in  direction  u  is  given  by 

(2)  W(K,u)  =  h ( K , u )  +  h (K , -u)  . 

A  convex  body  K  is  said  to  have  constant  width  b  if  and  only  if 
W(K,u)  =  b,  for  all  u  £  Sn_1 . 

For  a  plane  convex  body  K  we  shall  use  the  notation  h(K,0)  = 
h(K,u),  where  u  =  (cos  0,sin  0)  .  In  this  case  equation  (2)  can  be 
written  as 

(3)  W(K,0)  =  h (K , 0)  +  h (K , 0+*)  . 
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A  result  that  we  shall  find  useful  is  the  formula  of  Cauchy  for  the 
Euclidean  length  of  K,  namely 


n. 


From  (4)  we  can  obtain  Barbier’s  theorem  which  states  that  all  plane 
sets  of  constant  width  b  have  the  same  perimeter  »b .  An  ele;  entary 
proof  of  Barbier’s  theorem  is  given  in  Honsberger  [15] . 

The  area  of  K  is  denoted  by  A(K)  and  is  given  by 

r2v 

(5)  A  (K)  =  1  h(K,0)p(K,0)  60 

J0 

where  p(K,6)  denotes  the  radius  of  curvature  of  K  at  the  point  with 
outward  normal  u  =  (cos  0,sin  0)  .  As  a  consequence  of  the  Frenet- 
Serret  formulas  we  can  obtain  the  following  second  order  ordinary 
differential  equation  involving  radius  of  curvature  p(K,0)  and  the 
support  function: 

(6)  h  (K ,  0)  +  h(K,0)  =  p(K,0)  , 

where  dot  denotes  differentiation  with  respect  to  6.  Detailed 
discussion  of  (4),  (5)  and  (6)  can  be  found  in  Flanders  [12],  or  the 

monograph  on  convexity  by  Bonnesen-Fenche 1  [5] . 

For  a  plane  convex  body  of  constant  width  b  we  can  use  (3)  to 
obtain 

(7)  h (K , 0)  +  h(K,0+T)  =  b  . 

Substituting  0  by  6+n  in  (6),  we  obtain 

(8)  h(K,0+;r)  +  h(K,0+ff)  =  p(K,0+7r)  . 

Using  (7)  we  have 

3 


(9) 


h(K,0)  +  h  (K ,  0+*)  =  0  . 


Adding  both  sides  of  (6)  and  (8)  and  using  (9)  we  conclude  that  for  a 
plane  convex  body  of  constant  width  b  the  radius  of  curvature 
satisfies 

(10)  p( K,0)  +  p(K,0+v )  =  b  . 

Since  the  radius  of  curvature  p(K,0)  is  nonnegative  for  a  plane 
convex  curve,  we  use  (10)  to  obtain  (11)  for  a  plane  convex  body  of 
constant  width  b. 

(11)  0  <  p (K , 0)  <  b  . 

The  idea  of  using  optimal  control  is  simply  to  minimize  the  area 
A(K)  given  by  (5)  subject  to  differential  equation  (6)  and  conditions 
(7)  ,  (10)  and  (11)  . 

In  R3 ,  the  generalization  of  Barbier’s  theorem  is  the  fact  that 
the  total  mean  curvature  M(K)  for  a  set  of  constant  width  b  is 
constant.  That  is 

(12)  M  (K)  =  27rb  . 

We  also  have  the  following  remarkable  relationship  for  a  set  of 
constant  width  b: 

(13)  2V(K)  =  bS (K)  -  ^  b3  , 

where  V(K)  denotes  the  volume  of  K  and  S(K)  the  surface  area. 
Chakerian  [8]  contains  derivation  of  (12)  and  (13).  The  problem  of 
minimizing  the  volume  of  a  convex  body  of  constant  width  in  R3  is 
apparently  unsolved.  Using  (13),  minimizing  the  volume  is  equivalent 
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"to  minimizing  the  surface  area.  The  analogues  of  formulas  (5)  and 
(6)  for  volume  and  the  support  function  of  a  convex  body  in  R3  are 
given  in  Blaschke  [3] .  The  support  function  H  of  a  convex  body  in  R3 
sat i sf i es 

(14)  AH  +  2H  =  Rx  +  R^  , 

where  Rx  and  R2  are  principal  radii  of  curvature  and  A  is  essentially 
the  Laplacian  restricted  to  the  unit  sphere.  See  Bonnesen-Fenche 1 
[5]  . 

Chakerian  [9]  gives  the  following  lower  bound  for  n-d i mens i onal 
volume  V(K)  of  a  convex  body  K  of  constant  width  1  in  Rn: 


Chakerian  and  Groemer  [7]  give  an  interesting  survey  article  on 
sets  of  constant  width.  Bonnesen  and  Fenchel  [5] ,  Eggleston  [11] , 
and  Yaglom  and  Boltyanskii  [20]  have  good  treatment  of  required 
background  material  from  the  theory  of  convex  sets.  These  books  also 
contain  properties  of  sets  of  constant  width. 
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Optimal  Control 


In  the  -following,  we  sketch  the  theory  of  optimal  control  and 
give  Pontryagin’s  maximum  principle.  Boltyanskii  [4],  Leitmann  [18], 
and  Pontryagin  [19]  are  some  interesting  books  on  the  theory  of 
optimal  control.  Gelfand  and  Fomin  [13]  also  contains  an 
introduction  to  optimal  control .  Hermann  [14] ,  and  Brockett  [6] 
contain  differential  geometric  treatments  of  calculus  of  variations 
and  control  theory. 

We  will  consider  control  processes  that  can  be  described  by  a 
system  of  ordinary  differential  equations 

dx 

(17)  =  fj(x1,  .  .  .  ,xn,ux,  .  .  .  ,ur)  ,  l  =  1,2,.  ..,n 

where  Xj  ,  .  .  .  ,xn  are  space  coordinates  which  characterize  the  process 
and  u1,u2,...,ur  are  control  parameters  which  determine  the  process. 

In  order  to  determine  the  process  in  a  given  time  interval 
[to > tj] ,  it  is  sufficient  to  give  the  control  parameters  as  functions 
of  time  on  this  interval,  that  is 

|  (18)  uj  =  Uj(t)  ,  j  =  l,...,r  . 

Assuming  that  the  problem  is  well  posed,  for  a  given  initial 
s  te 

(19)  xj(t0)  =  xj°  , 


the  solution 
Consider 


(20) 


of  (17)  is  uniquely  determined, 
the  functional 

J  —  I  fo(x^,...,xn,U2,..-,Ur)  dt  ■ 
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For  each  control  (18)  on  [tojtjJ  ,  the  process  is  determined  and  the 
■functional  J  assumes  a  certain  value.  Assuming  that  there  is  a 
control  (18)  which  transfers  the  object  from  a  given  initial  state 
(19)  to  a  final  state 

(21)  Xj(ti)  =  xj1  , 

the  object  is  to  find 

(22)  uj(t)  ,  j  =  1  ,  .  .  .  ,  r 

which  transfers  the  object  from  (19)  to  (21)  in  such  a  way  that  the 
functional  (20)  has  a  minimum.  In  general  there  are  restrictions  on 
control  parameters  Uj .  Thus  we  shall  assume  the  vector  u  belongs  to 
a  region  U  in  2-d imens i onal  Euclidean  space  called  the  control 
region ,  that  is 

(23)  u  €  U  . 


An  admissible  control  is  an  arbitrary  piecewise  continuous 
control  in  the  control  region  U. 

We  now  state  Pontryagin ’s  maximum  principle.  As  a  good  start  for 
understanding,  one  may  study  the  shortest  time  problem  for  a  phase 
point  moving  from  an  initial  point  x°  to  the  origin  in  accordance 
with  =  x2 ,  =  u,  ju|  <  1.  This  problem  is  studied  in  many 

books  including  Boltvanskii  [4] . 


Maximum  Principle 

Let  u(t)  ,  tG  <  t  <  tj ,  be  an  admissible  control  such  that  the 
corresponding  trajectory  x(t)  which  begins  at  the  point  x°  at  the 
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time  tQ  passes  at  the  time  t,  through  a  point  x1 .  In  order  that  u(t) 
and  x(t)  be  optimal  it  is  necessary  that  there  exists  a  nonzero 
continuous  vector  function  ^(t)  =  (^0 (t)  ,  .  .  .  , (t ) )  corresponding  to 
u(t)  and  x(t)  satisfying 


(24) 


d^i  _  _  dH 
cTt  dxj 


i  =  0,1,2,. ..,n, 


where 


(25)  H(^,x,u)=£  ipQfa  , 

a=0 

such  that 

(a)  for  every  t,  t0  <  tj  ,  the  function  H  ( tp  (t )  ,  x  ( t )  ,  u )  of  the 
variable  u  €  U  attains  its  maximum  at  the  ponit  u  =  u(t) : 

(26)  H(^(t)  ,x(t)  ,u(t)  )  =  m(v!>(t)  ,x(t)  )  =  sup  H(^,x,u)  . 

ueu 


(h)  At  the  terminal  time  tj  ,  the  relations 
(27)  iZ’o(ti)  <  0,  m(t/>(t1)  ,x(tj) )  =  0 


are  satisfied.  It  turns  out  that  if  V’(t),  x(t)  and  u(t)  satisfy 
systems 


(28) 


^Li  _  dH 

dt  <9xj 


i  =  0,1, . . . , n  , 


(29) 


M  _  _  M 

dt  dXj  ’ 


i  =  0,1,.. . , n  , 


and  condition  (a),  then  the  time  functions  </>Q(t)  and  m(V’(t )  ,  x  (t)  )  are 
constant.  Thus  (27)  may  be  verified  at  any  time  t,  t0  <  t  <  t1?  and 
not  just  tj  . 
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Formu  1  at:  ion 


The  Blaschke-Lebesgue  theorem  states  that  the  Reuleaux  triangle  has 
the  least  area  of  all  plane  convex  sets  K  of  the  same  constnat  width 
b.  The  minimum  area  is  - T  .  In  the  following,  without  loss  of 

generality,  we  assume  b  =  1  .  We  also  write  our  formulas  concerning  a 
plane  convex  body  K  without  specifying  K.  For  example  (6)  will  be 
written  as 


(30)  h(0)  +  h  (0)  =  p(0)  . 

Thus  we  can  rewrite  formulas  (5),  (6),  (7),  (10)  and  (11)  as  follows: 


(31) 

A  = 

1 

2 

h(0)p(0)  d6 

Jo 

(32) 

h(*) 

+ 

h(0)  =  p(0 )  , 

(33) 

h(0) 

+ 

h(0+^)  =  1  , 

(34) 

P(0) 

+ 

P(0+ *)  =  1  , 

(35) 

0  < 

p(0)  <  1  • 

We  also  have  the  following  formula  for  any  convex  body  K  which 
can  be  derived  from  Frenet-Se r ret  formulas  (Flanders  [12]). 

h(0)  =  x(0)  •  t(0), 

where  x  denotes  the  position  vector  and  t  unit  tangent  vector  at  the 
point  where  the  outward  unit  normal  is  given  by  (cos  0, sin  0) . 

Through  the  endpoints  of  any  diameter  of  a  set  K  of  constant 
width  there  are  support  lines  of  K  perpendicular  to  that  diameter. 
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Furthermore  there  is  a  diameter  for  which  the  corresponding  support  j 
lines  are  tangent  to  the  curve  (Eggleston  [11] ,  P.  126).  Let  one  ' 
endpoint  of  such  a  diameter  be  taken  as  origin.  Using  (33)  and  h (0)  j 
we  obtain 


(36)  h (0)  =  1,  h(T)  =  0  , 

(37)  h(0)  =  0  ,  h(T)  =  0  . 
Substitute  (33)  and  (34)  in  (31)  to  derive 


(38) 

A  =  l  f  [1  +  2h  (0)p(0)  -  h(0)  -  p(0)] 
J0 

d&  . 

Let 

(39) 

xl(*)  =  2h  (0)  -  1  , 

(40) 

x2(0)  =  xx(0)  =  2h(0),  and 

(41) 

u  (9)  =  2p(0)  -  1  . 

If  we  now  substitute  (39),  (40)  and  (41)  in  (32), 

(35),  (36).  (37) 

and  (38) 

,  we  obtain  the  following  approach  to  proving  the  Blaschke- 

Lebesgue 

theorem : 

(42) 

fir 

Minimize  ^  (1  +  x1(0)u(ff))  d0  , 

J0 

subject 

to 

(43) 

X  j  —  X2  » 

(44) 

X  2  =  U  —  X-1  , 

(45) 

xl(0)  =  1,  x!  (ir)  =  -1  , 

10 


(46) 


x2(0)  =  0,  x2(ir)  =  0  , 


(47)  | u |  <1  . 


Analysi s 


The  Blaschke  selection  theorem  states  that  every  infinite 
sequence  of  closed  convex  subsets  of  a  bounded  portion  of  Rn  contains 
an  infinite  subsequence  that  converges  to  a  closed  nonempty  subset 
(Eggleston  [11],  P.  64).  The  Blaschke  selection  theorem  implies  that 
the  minimum  area  exists. 

Minimizing  (42)  subject  to  (43)-(47)  is  equivalent  to  minimizing 


(48) 


xx(0)u(0)  dO 


under  the  same  constraints.  We  proceed  to  use  Pontryagin’s  maximum 
principle  to  minimize  (48). 

Let 


(49)  ^  =  XjC^uCtf)  . 

Let  =  (^o(^)  >V’i(^)  i^2(^))  be  the  auxiliary  vector.  Use  (25)  to 

obtain 


(50) 


H  =  t/>0X1U  +  ^2x2  +  ^'2  (u-xl)  =  (V’oxl+V,2) u  +  Vtx2  -  */>2xi 


We  then  use  (24)  to  write  down  differential  equations  for  i> . 


(51) 


dip0  _ 
dO  ~ 


0H  _  n 

ET0  ~  0  ’ 


(52) 


dF  -  '  ax,  ~  ’ 


l 


l 


ii 


d^2  _  3H  _  . 

so  -  ~  &r2  - 


There  exists  a  nonzero  continuous  vector  ip(0)  satisfying  (51)-(53) 
such  that 

(54)  max  H  (x (0)  ,  u  ,ip(0)  )  =  H  (x  (0)  ,  u  (0)  ,  \fc(0)  )  . 

M<1 


Now  we  consider  cases  to  analyze  the  maximum  of  H  given  in  (50) 
as  a  linear  function  of  u. 

Case  1(a):  rp0  0  and  ip0*-i  +  ip2  ^  0  (not  identically  equal  to 
zero)  lead  to  a  contradiction. 

In  this  case 


f  1  if  V’oXi  +  ^2  >  0  , 

u  =  { 

(-1  if  iPo*-i  +  ^ 2  <  0  . 


Using  (43)  and  (44)  we  obtain  xx  =  u  +  A  sin(0-a)  for  some  A  and  a. 
Continuity  of  ipQ ,  Xj  and  i<2  imply  that  ^0xx  +  \p2  =  ^o(u  +  A  sin(r-o)) 
+  ip2  is  continuous  at  a  switching  point  r.  Continuity  at  r  implies 

(56)  t/>0  +  t/’o  A  sin(r-a)  +  ip2(r)  =  -^o  +  V'o  A  sin(r-a)  +  ^’2(0 


Hence  using  (56)  we  conclude  that  ip0  —  0  which  is  a  contradiction 
Case  1(b)  :  Assume  ip0  £  0  and  0oxi  +  t/>2  =  0  • 

In  this  case  tp2  =  -ip oxi .  Hence  ip2  =  -ip0*-i  =  ~*l>i  =  ox2  >  and 
(57)  =  ip2  -  ip0u  =  iPo'X-2  =  ipo^u-yi^  . 


Hence 


-^oxi  -  V'ou  =  V'oU  -  t^oxi 
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which  implies  u  =  0.  Using  (41)  we  conclude  p(0)  =  ^  which 
corresponds  to  a  circle  of  radius  giving  the  maximum  area  rather 
than  the  minimum. 

Case  2 .  Suppose  tp0  =  0 . 

In  this  case  ip2  ^  0.  Since  ip2  =  0  and  (53)  imply  Vj  =  0.  But 
we  know  that  V  is  a  nonzero  vector.  Hence  using  (55)  we  obtain 


(58) 


if 

if 


Differential  equations  (52) 


^2  >  0 
ip  2  <  0 

and  (53) 


♦ 


reduce  to 


(59) 

(60) 


dip  i  . 

W  =  *2 


dt/>2 


d0 


=  -V'i 


An  analysis  similar  to  Pontryagin  [19,  pp .  27-35]  will  give  the 
switching  curve  in  Figure  2  where 


(61) 


above  the  cuve 

below  or  on  the  curve  . 


We  can  now  interpret  the  system  of  differential  equations  (43) 
and  (44)  as  equations  of  motion.  Our  objective  would  be  to  get  from 
(1,0)  to  (-1,0)  in  such  a  way  that  we  minimize  (48).  One  can  show 
that  this  is  equivalent  to  minimizing  the  action  of  the  moving 
object . 

Optimal  trajectory  starts  at  (1,0).  Since  (1,0)  is  above  the 
curve  in  Figure  2,  we  use  u  =  -1  until  the  trajectory  intersects  the 
switching  curve  on  the  x2  axis.  We  will  then  use  u  =  +1  which  will 
lead  to  (-1,0).  If  we  interpret  0  as  time,  then  the  total  time  from 
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- 1  are 


(1,0)  to  (-1,0)  is  We  realize  that  xx  =  1  and  xx  = 

solutions  of  (43) -(46)  in  neighborhoods  of  9  —  0  and  9  =  r 
respectively.  Hence  in  order  to  use  the  total  time  suggested  by 
functional  (48),  the  object  starts  at  (1,0)  and  waits  there  for  time 
a  <  Then  the  object  goes  to  (0,-^3)  and  then  to  (-1,0).  During 

the  time  interval  (^  +  a,*)  the  object  waits  at  (-1,0)  for  a  period 
of  g-a.  In  this  way  the  total  time  would  be  ir .  See  Figure  3. 

We  use  differential  equations  (43)-(46)  and  (62)  below  to  obtain 
(63). 


(62) 


u  = 


1 , 

0  <  9  < 

a 

-1 , 

a  <  9  < 

§  +  a 

1 , 

3  +  *  <  * 

<  ^  +  a 

-1 , 

a  +  ^  < 

VI 

(63) 


h(»  =  1  tail) . 


1  ,  0  <  9  <  a 

cos  (9 -a)  ,  a  <  9  <  g  +  a 

1  +  cos(0  +  ^-a),  5  +  q  <  0  <  ^  +  a 

0, 


s 


<  IT 


However  (63)  gives  the  support  function 
width  1 .  The  angle  a  corresponds  to  the 
Reuleaux  triangle  by  angle  a  will  result 
where  6  is  replaced  by  a  +  g,  0  <  a  <  £. 
obtain 


of  a  Reuleaux  triangle  of 
fact  that  a  rotation  of  a 
in  a  new  support  function 
We  can  now  calculate  (42)  to 
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(64) 


i  f2ir 

Minimum  area  =  4  (1  +  x1(tf)u(^))  dO  = 

J0 


_  *  -  ^3 
2 


Hence  by  comparison  of  case  1  and  case  2  we  choose  the  minimum,  which 


is 


as  desired. 
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Figure  2.  Switching  Curve 
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